arXiv:1505.05307v2 [nucl-th] 21 Jan 2016 


Effects of chiral imbalance and magnetic field on pion superfluidity and color superconductivity 


Gaoqing Cao and Pengfei Zhuang 

Department of Physics, Tsinghua University and Collaborative Innovation Center of Quantum Matter, Beijing 100084, China 

(Dated: January 22, 2016) 

The effects of chiral imbalance and external magnetic field on pion superfluidity and color superconductivity 
are investigated in extended Nambu-Jona-Lasinio models. We take Schwinger approach to treat the interac¬ 
tion between charged pion condensate and magnetic field at finite isospin density and include simultaneously 
the chiral imbalance and magnetic field at finite baryon density. For the superfluidity, the chiral imbalance 
and magnetic field lead to catalysis and inverse catalysis effects, respectively. For the superconductivity, the 
chiral imbalance enhances the critical baryon density, and the magnetic field results in a de Haas-van Alphan 
oscillation on the phase transition line. 
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I. INTRODUCTION 

Quantum systems in an external magnetic field exhibit var¬ 
ious unexpected and interesting features. In condensed matter 
physics, quantum hall effect and fractional quantum hall ef¬ 
fect were discovered in two-dimensional electron systems in 
1970s and 1980s UHl. More recently, quantum spin hall ef¬ 
fect was proposed to exist in special two-dimensional electron 
systems with strong spin-orbit couplings ail, and quantum 
anomalous hall effect was found in a ferromagnetic material 
chromium-doped (Bi,Sb)2Te3 in 2013 H- For a dense non- 
relativistic system, the de Haas-van Alphan oscillation was 
observed in some thermodynamic systems ai- In high en¬ 
ergy physics, this oscillation is found in nuclear and quark 
matters iMl- It is beyond our expectation that the mag¬ 
netic field reduces the critical temperature of chiral symmetry 
restoration in Quantum Chromod yna mics (QCD), known as 
the inverse magnetic catalysis |[l3l - [0 . 

In high energy heavy ion collisions, there might exist 
charge asymmetry on the opposite sides of the reaction plane, 
induced by the chiral magnetic effect ifisl - l^ . Besides a very 
strong magnetic field (up to eB ~ IQml at the LHC energy) 
created in the early stage of heavy ion collisions, the randomly 
generated chiral imbalance plays also an important role in the 
hot medium. In neutron stars where the magnetic field and 
baryon chemical potential are both very large, chiral current 
can be induced along the magnetic field by triangle anoma¬ 
lies ll23l - l2^ . The chiral current will eventually cause chiral 
separation and generate chiral imbalance which can be quite 
large in the two hemispheres along the magnetic field. While 
the chiral imbalance or instead chiral chemical potential 
together with the magnetic field B created by the spectators 
have been widely investigated in the study of chiral symme¬ 
try restoration 11261 - 1^ at finite temperature, it is still an open 
question how the chiral imbalance and magnetic field affect 
pion superfluidity at finite isospin density and color supercon¬ 
ductivity at finite baryon density. 

The pion superfluidity in an external magnetic field was 
studied in Nambu-Jona-Lasinio (NJL) model S, linear 
sigma model and Lattice QCD 1^ . When charged pion 
condensate exists, u and d quark fields are no longer the eigen¬ 
states of the quark propagator in flavor space, and a direct 
extension from the case without magnetic field B to the case 


with finite B by substituting the transverse momenta with the 
Landau levels of u and d quarks is invalid. The model cal¬ 
culation at hadron level shows a magnetic catalysis which is 
inconsistent with the lattice simulation at quark level. The 
color superconductivity in an external magnetic field is not so 
complicated as pion superfluidity. While the color conden¬ 
sates are also charged, the effective charges of u and d quarks 
with red (green) and green (red) colors are opposite to each 
other d, when one takes the interaction between quarks and 
massless gluons into account. Thus, the color condensates can 
be effectively treated as neutral in an external magnetic field. 

We investigate in this paper both the chiral imbalance and 
magnetic field effects on pion superfluidity and color super¬ 
conductivity in extended NJL models. Considering the diffi¬ 
culty to treat the charged pion condensate in a magnetic field, 
we study separately the chiral imbalance effect and magnetic 
field effect, and for the latter we employ the Schwinger ap¬ 
proach to include the interaction between the charged pion 
condensate and the magnetic field. For the color superconduc¬ 
tivity, we take into account the chiral imbalance and magnetic 
field simultaneously. 


II. PION SUPERFLUIDITY WITH CHIRAL IMBALANCE 


The Lagrangian density of the extended NJL model with 
quark isospin chemical potential pi and chiral chemical po¬ 
tential ps is defined as 

£. = ij/[i$ - mo ■+ yroTs + P57075j^lt+G [{ifnlff + {^i75Tt//fj , 

( 1 ) 

where \f/ - {u, dY is the two-flavor quark field, t , are pauli 
matrices in flavor space, otq is the current quark mass, and 
G is the coupling constant with dimension GeV“^. In order 
to study the bound states of the system, we introduce four 
auxiliary fields cr = -2G4nj/ and ti - -2Gi[/i7sTtl/, and the 
Lagrangian density becomes 


+B57075 


Pi 

id-mo-cr- ijs {t^tio + T±n±} + —70T3 


-A- 


(Y + 7:^+ n^n± 
4G ’ 


( 2 ) 
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where the auxiliary fields are related to the physical fields cr, 
ttq = 713 and Ti± - (tti + ini) / V2, and t± = (ti + ixa) / V2 are 
the raising and lowering operators in flavor space. 

The order parameters of spontaneous chiral symmetry 
breaking and isospin symmetry breaking are respectively the 
expected values of the auxiliary fields (cr) and {n±) ((no) - 0). 
Taking them to be real and constants, considering the relation 
between the chiral order parameter and the dynamical quark 
mass (cr) = m — nio, and denoting (7r±) by A, the thermody¬ 
namic potential of the quark system in mean field approxima¬ 
tion can be expressed in Euclidean space as 


chemical potential, the pion superfluidity phase with A 9^ 0 is 
always the ground state of the system 11^ . 

In vacuum, the critical coupling constant to keep sponta¬ 
neous chiral symmetry breaking is Gc - n^l{2NcA^). By 
scaling the coupling constant G by Gc, the gap equation for 
the nonzero pion condensate is reduced to 



where p,iii,^s and A are all scaled by the cutoff A. 


Q. = 


(m - mof + A^ 


4G 


(3) 


^Trln - m - /ysTjA -i- yroTs + /tsyoTs) 


with = iyoPo - JiPi and f) - l/T, where the trace is taken 
over the quark spin, flavor, color, coordinate and momentum. 
At zero temperature, the potential is explicitly expressed as 


Q = 


{m — mo)^ + A^ 
4G 



(4) 



where Ei j(p) - yj^fjip) + ^ are the quasiparticle energies 

with ^ijip) - yjip + ips)^ + np- + jpijl, and Ac = 3 is the 
color degrees of freedom. The two order parameters cr or m 
and A are determined by the minimum of the thermodynamic 
potential, dClIdm - 0 and dQldA = 0, which lead to the two 
coupled gap equations. 


{m - mp) 
2G 
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/V y f _ m _ 

!j=± J VdPi + ipiY + 


f 
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A 
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(5) 


We now discuss the effect of chiral imbalance on chiral 
symmetry restoration and pion superfluidity. Since the NJL 
model is not renormalizable, it is necessary to introduce a mo¬ 
mentum cutoff A to regulate the integrations in the gap equa¬ 
tions 0. Let us first consider the possible phases of the sys¬ 
tem in chiral limit with vanishing current quark mass mo = 0. 
In this case, there are only two parameters G and A in the 
model, they are fixed to be G = 5.01 GeV“^ and A = 0.65 
GeV by fitting the chiral condensate and pion decay constant 
in vacuum with T - pi — ps — ^ ll^ . As can be seen from the 
gap equations 0, there always exists a trivial solution m - 0 
and A = 0, corresponding to the phase with chiral symmetry 
and without pion condensation. The nontrivial solutions are 
usually preferred by the system, as they correspond to lower 
thermodynamic potentials. However, the two gap equations 
with m, A 0 contradict with each other ll^ . and there are 
only two possible nontrivial solutions; the phase with spon¬ 
taneous chiral symmetry breaking characterized by m 0 
and A = 0 and the phase with spontaneous isospin symme¬ 
try breaking described by A 0 and m - 0. At finite isospin 


FIG. 1: (color online) The pion condensate A as a function of isospin 
chemical potential p; at fixed chiral chemical potential p 5 in chiral 
limit. All the quantities are scaled by the momentum cutoff A. 



FIG. 2: The pion condensate A as a function of chiral chemical po¬ 
tential p 5 at fixed isospin chemical potential p/ in chiral limit. All the 
quantities are scaled by the momentum cutoff A. 

The scaled pion condensate is shown in FiglT]as a function 
of scaled isospin chemical potential. As a first order phase 
transition, the pion condensate A jumps up from zero to a fi¬ 
nite value at the critical isospin chemical potential p) = 0. 
While the condensate increases with increasing p/ or ps in 
general case, it drops down with increasing ps when p/ is large 
enough, see the up-right corner of FiglT] Considering the fact 
that Pi in this case is already close to or even beyond the cut¬ 
off A, this dropping down is probably an artifact of the model. 
The p 5 dependence of the condensate at fixed p//A = 0.4 is 
shown in Fig|3 which displays a monotonous increase. Note 
that, while the pion condensate increases with chiral imbal¬ 
ance at reasonable isospin density, the critical point of pion 
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superfluid is not affected by the chiral imbalance, it is always 
located at /i/ = 0 . 



FIG. 3: (color online) The chiral and pion condensates m and A 
as functions of isospin chemical potential p; in real case. All the 
quantities are scaled by the momentum cutoff A, and the solid, 
dashed and dotted lines correspond to fixed chiral chemical poten¬ 
tial ps/A = 0, 0.2 and 0.3. 



yUs/A 


FIG. 4: The chiral and pion condensates m and A as functions of 
chiral chemical potential /J 5 at fixed isospin chemical potential /j; in 
real case. All the quantities are scaled by the momentum cutoff A. 

In real case with finite current quark mass, the three param¬ 
eters of the model are fixed to be G = 4.93 GeV“^, A = 0.653 
GeV and niQ - 5 MeV by fitting the pion mass m„ - 134 
MeV, pion decay constant G = 93 MeV and quark condensate 
(cr) - -2 X (0.25GeV)^ 1^ . The isospin chemical potential 
dependence of the two order parameters is shown in Figl3] At 
fis/A = 0, the superfluid starts at the critical isospin chemical 
potential= m„ (/rJ/A = 0.21) 1^ . With increasing yus, the 
critical point shifts towards the left but the pion condensate 
increases, which indicate a catalysis effect. The quark mass 
keeps as a constant in the normal phase at /j.j < //J and drops 
down monotonously in the pion superfluid. With increasing 


P 5 , the quark mass increases in the normal phase but decreases 
in the superfluid. Suppose the chiral phase transition happens 
at the same critical point as the pion superfluidity, the /us de¬ 
pendence of the quark mass shows an inverse catalysis effect 
on the chiral phase transition. This is similar to the lattice sim¬ 
ulated magnetic field effect at finite temperature CMl: The 
chiral condensate increases at low temperature but the critical 
temperature is reduced. The chiral imbalance effect on the two 
order parameters at fixed isospin chemical potential is clearly 
shown in FigH] 


III. PION SUPERFLUIDITY WITH MAGNETIC FIELD 


Now we turn to the pion superfluidity in an external mag¬ 
netic field. The Lagrangian density of the NIL model is writ¬ 
ten as 

Z = if/ [iU) - mo + y70^3) ^ + G ] , ( 7 ) 


where = dfj+iqA^^ is the covariant derivative in flavor space 
with electric charges qu - 2e/3 and qa - -e/3 for u and d 
quarks, and the potential - ( 0 , 0 , Bxi, 0 ) defines a constant 
magnetic field along the X 3 -axis through B = V x A. Since 
it is cumbersome to study pion superfluidity in an external 
magnetic field in quark models, as mentioned in the introduc¬ 
tion, we derive here the phase transition line of pion superfluid 
by using the Ginzburg-Landau (GL) approach. Following the 
same procedure as in SectionUT] the thermodynamic potential 
in Minkowski space is expressed as 


Q = 


(m — nio)^ + A^ 
4G 


-I- —Trln 
V 4 


(iGu) ' -iysA \ 
-iysA* (/Grf)-'j’ 


( 8 ) 


where V 4 is the space-time volume, and {iGf)~^ = i0f - m± 
(juil2)yo if = M, d) is the inverse quark propagator at mean 
field level. We now take Taylor expansion of Q in terms of the 
pion condensate A around the critical point and keep only the 
first two terms. 


^ 4y^TrG„iGyi3'2^.v2y3^3' 

= JIA^ + 

2 


(9) 


with the quark propagator in coordinate space 


Gxy - 


I iGu{x,y) 

0 \ 

1 ® 

-175 Ay \ 

1 0 

iGd(x,y) j 

\ - 175 A; 

0 I 


( 10 ) 


and the coefficients ^ and S in the GL approximation 
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^ ^ ^ + ■^Tr[G„(x,y)!y5G^(y,x)/75e “'I., 

S = --^T'r[G,Xx,yi)iy5Gd{yi,y2)i75Gu{y2,y3)i75Gd(y3,x)iy5e~“'i^^^- (n) 


Note that in the expression of the coefficients, we have taken into account the interaction between the charged pion condensate 
and the magnetic field which is not included in the NJL model Q- It is introduced by a straight-line link (the exponential in 
(fTTb ) between two points Ay and A*, and each condensate is linked only once. The quark propagators Gf(x,y) can be evaluated 
with Schwinger approach 11^ . 


Gf(x,y) 

Sf(x) 


!y 

-■/ 


Sf(x- 
ds 


■y\ 


( 12 ) 


I6(7rsy 


e '[^'”"+4,(4 4 [cotB^ - 7172] m H- ^ (/q - ^3 “ B } ((^1 -H ^2) cotB;^ H- ^21 - ^12)) 


with Bj = q/sB, = j^Xy and x]_ = x^ + x^, where A^. = {+fj.jl{2qf), 0,0,0) + A^ are effective potentials for u 

and d quarks, and the integration in the exponential from 7 to x is along a straight line. With the known quark propagators, the 
coefficients can be expressed in terms of the functions S /, 


^ ^ J d\x - 7)Tr [e'T(-o-.vo) 5 „(x - y)ir 5 e-‘"^^^«~^°^Sd(y - x)ij 5 ], 

S = r - 72 )t /'^(72 -73)x (13) 

^4 J 

Tr l^e'f _ y^)iJ^e^‘ib•lo-yw)g _ y2)f75e'^*^-““^™^51,(72 - y 3 )ij 5 e~‘’^^^'’~'‘''^S 1/(73 - x)f75], 


where the Wilson lines in the quark propagators are exactly canceled by the interaction between the charged pion condensate 
and the magnetic field, and therefore the gauge invariance is guaranteed. For convenience we transfer from the coordinate space 
to the energy-momentum space and switch on the temperature. The quark propagators in Euclidean space are then given as ll^ 


poo ^ 

Sf(a>n,k) = -i I ^ _ A: 2 i) tanh) (l - / 7 i 72 tanhB}) 


(14) 


with the Matsubara frequency = (2n + \)nT (n 6 Z) of fermions. From the GL theory, the condition ^ - 0 determines the 
second order phase transition line of the pion superfluidity. After a straightforward calculation, we obtain 


1 r c/^k 

^ =- NcT > -^Tr 

4G " 4^J (2;r)3 


Su{(^n+ ’ kj iysS^ [oJn - *y - kj iys 


(15) 


^ 4A,7’ ^ J ^ JdsJ dfR(s, t, oj„, k) (1 + fMMO) + ki (l - f^(s)) (l - 


4G 


with fu(s) - tanhB*, fj(t) = tanhB^, transverse momentum k 5 _ = kj H- and 


"Ris t to,, k) = 


(16) 


where the integrations over s, t and k are divergent and a regularization scheme is needed. We take the way given in Refs, ll^l^ 
and absorb all the divergence into the vacuum term which is then regularized by the three-momentum cutoff A. In this way the 
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convergent coefficient ^ includes a vacuum part and two magnetic field dependent parts J{.bi and 


- — + .^0 +-^Bl +-^B 2 , 

40- 


-^0 - 




7.71^ 
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/Bf.+fd{ 


IB', 


(17) 


with EiC't) = + /r/^ + //i//2. ft should be remarked that there is an implicit condition m > jiijl here, which makes the 

integrations over s convergent in ultraviolet domain and which is guaranteed at small pion condensate where the dynamical 
quark mass m is relatively large. This condition is only due to Schwinger approach itself for introducing proper-time integral 
and there is no such restriction with Ritus method 

The derivation of the coefhcient S is more tedious. After a careful calculation, we finally obtain its explicit expression. 


S = ANcT 


( 18 ) 


y X (2 V X , t', (On, k) X 

1 [(u 7 >f + k]+ + (mfxif + (k3/2/)2] y (1 + //„(.)) (1 + //„(/)) (1 + IMt)) (1 + IMt')) 

+y (l - /«(^)) (l - (l - /1(0) (l - + ki [4 (cu>^ + kl + m^)il + fu{s)m) (l - /;?(s')) (l - fAt')) 

- [{c 4 f + kl + ,A) (1 - fu{s)fu{s')) (l - fAt)) (l - fAt')) - (^{cA,f + kl + 2772) (1 - MDfdit')) (1 - / 2 (s)) (1 - fAs')) \ 


with the effective frequencies 02 “ = a;,, + // 2//2 and o/f = w„ - /A 2 // 2 . ft is not necessary to numerically calculate the coefhcient S, 
if we only focus on the phase transition line of the pion superfluid. What we are interested in here is its sign which determines 
the order of the phase transition and thus the validity of the GL approach. Keeping in mind the symmetry between s{f) and s'{f) 
in the integrations, we have the inequality 

4(1+ mfM) (1 - fA^')) (1 - fAt')) - (1 - fuis)fu{s')) (1 - fAt))[i - fAt')) - (1 - fMAit')) (i - fA^')) (i - f{s)) 

> (1 - fui^')) (1 - fAt')) [4(1+ fu{s)fm - (1 + fA^)) (1 - fAt)) - (1 + fAt)) (1 - f(^))] 

= 2(1+ fu{s)fAt)f (1 - y(*')) (1 - fAt')) > 0. (19) 

Then using the condition m > fj.jl2 for small condensate A and the fact of small ratio Im K/Re 7?, the coefficient S is found to 
be positive dehnite around the transition line of the pion superfluid. Thus, the transition from normal phase to pion superfluidity 
with non-vanishing magnetic field is proved to be of second order. 

The phase transition line of the pion superfluid ^ - 0 is coupled with the dynamical quark mass m which should be evalu¬ 
ated consistently with the gap equation for the chiral phase transition. Following Schwinger approach and taking the vacuum 
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regularization scheme adopted above, we have the gap equation 

0 = 


^ f=u,d 


2G 

m — niQ Ncm? 


2G 

Nctn 

Ati^ 


71 “ 

ds 


A 2 


A -v/1 H- - - min 
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A / A2 

- + a /1 + — 




Si 


Iddh 


V^2 + ot2 1 + e^Ak)/T 


r” ^ ./ 2 / ^_L_\(q„sB -Bi — 4 \ 


T^\/^dsB 

Ms) 


- 1 


( 20 ) 


where d^(z,q) is the third Jacobi theta function obtained by on the phase transition line. By comparing it with the upper 
working out the summation over the Matsubara frequency. panel, there is always m^ > /rj/2. 
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FIG. 5: (color online) The phase diagrams of pion superfluid in p; - B 
plane (upper panel) and the quark mass on the phase transition 
line (lower panel) in chiral limit (mo = 0) and real case (mo = 5 
MeV). 

We now show the phase diagram of pion superfluidity in an 
external magnetic field. We adopt the same set of parameters 
of the model as given in Section|IIl The phase diagram in p/ - 
B plane at zero temperature is shown in the upper panel of 
FiglS] The two solid lines are the phase transition lines of pion 
superfluidity in chiral limit with mo = 0 and real case with 
mo = 5 MeV. The normal phase without pion condensate is 
under the corresponding line, and the pion superfluidity phase 
is above the line. In both cases, the critical isospin chemical 
potential ffj increases with the strength of the magnetic field B, 
which indicates clearly the effect of inverse magnetic catalysis 
on the pion simerfluidity. This is consistent with the lattice 
QCD result IBlI] . Note that the above analytic and numerical 
results depend on the condition m > p//2. To check if this 
condition is satisfied, we show in the lower panel of Figl^the 
dynamic quark mass nic as a function of the magnetic field 


FIG. 6; (color online) The phase diagrams of pion superfluid in p/ - T 
plane (upper panel) and the quark mass m^ on the phase transition 
line (lower panel) at fixed magnetic field and in real case with mo=5 
MeV. 

The phase diagram in p/ - T plane is shown in the upper 
plane of Fig|6]in real case with mo=5 MeV. The phase transi¬ 
tion line at a fixed magnetic field separates the normal phase 
at low p/ and the pion superfluidity at high p/. The increas¬ 
ing critical isospin chemical potential pj with the magnetic 
field indicates again the inverse catalysis effect on the pion 
superfluidity. While the temperature dependence of p)^ is very 
smooth at low temperature, it goes up quickly at high temper¬ 
ature. The quark mass on the phase transition line is shown 
in the lower panel of Fig|6] It decreases with temperature due 
to the gradual chiral symmetry restoration. The small trian¬ 
gle structure in the intersection region shows the de Flaas-van 
Alphan oscillation 10,0] induced by the interplay among the 
field, chemical potential and temperature. The reason why 
the de Haas-van Alphan oscillation doesn’t show up at low 
temperature is that the gap equation (l20l i only depends on p/ 
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weakly now as the third Jacobi theta function 1 in 

the forth term on the right hand side. 


IV. COLOR SUPERCONDUCTIVITY WITH CHIRAL 
IMBALANCE AND MAGNETIC FIELD 


We discuss now color superconductivity including simulta¬ 
neously the chiral imbalance and magnetic field effects. The 
Lagrangian density of the extended NIL model with baryon 
chemical potential fig, chiral chemical potential //5 and exter¬ 
nal magnetic field B is defined as 

-C ^ ifr [iy^ [d^ - igT^Gl) - mo + jisyo + Bsyoys] ( 21 ) 

+Gs + Gd (iiffcseoystf/) {iij/eeiysipc), 


partition function 


Z = 


J" [dif/] exp ^ - i J'd' 


4Gs 

■f, - m + jj.Byo + Bsyoys) ^ 


{m — mor 


4Go 


+ -^c (iiffsesysil/c + iil/cseoystf/) 


(25) 


where the functional integrations over ^ and can be accom¬ 
plished in the Nambu-Gorkov space. With the help of the 
charge projectors 




' diagiO,0,0,1,1,0), 
diag(0,0,0,0,0,1 ), 
^f/ag(l, 1 , 0 , 0 , 0 , 0 ), 
diag(0,0,1 ,0,0,0), 


q^-l/2 

q^O 

q^Ml 

q^\. 


(26) 


where i/^c = and tf/c - tf/^C are charge-conjugate spinors 
with C = iyiyo, Sij and {e 3 )ab = are respectively the an¬ 
tisymmetric matrices in flavor and color spaces with indexes 
i,j = {u,d) and a,b = {r,g,b), the term -yigAsG^tf/ll ac¬ 
counts for the interaction between quarks and massless glu¬ 
ons with Tg being the 8th Gell-Mann matrix, and 

the coupling constants Gs and Go in the scalar and diquark 
channels are related to each other by the Fierz transformation 
Gd = 3Gs/4 m. 

Using the ’’rotated” charge operator Q - Q ® Ic - 1/ ® 
/lg/(2 V3), a massless Gem(l) held = A^cos0 + G^sin0 
and a massive gluon field G* = G^ cos 0 - A^ sin 6 are ob¬ 
tained d. By neglecting the massive gluon field and re¬ 
placing the massless held with an external magnetic held 
A^ = (0,0, Bxi, 0) as 0 is small, the following simplified La¬ 
grangian is derived, 

-C - ^ [iy^- mo + lUByo + jusyors] ^ (22) 

+Gs + {ii/iysT^f\ + Gd {i^cse^ysy) (iiffseoystf^c) 


the bispinors in Nambu-Gorkov space are defined as 


= (27) 


with - Pqifi , and the partition function becomes 


= J" [dtp] [dip] exp I - / J" 


(m — mo)^ 


4Gs 


4Gd 


(28) 


where the inverse fermion propagators in Nambu-Gorkov 
space are defined as 


Sq^ 


i(D;)-' 0 \ 

0 /(Dp-I I 

-iAcT2A2y5P-q \ 
-iAcT2A2y5Pq i(Dq)~^ ) 


q^0,l 

(29) 


with Df,^df, + iQeAf, and Q = diag(q^,, q^g, q^o, q±-, qdg, qdb) 
- diag{l/2,l/2,l,-l/2,-l/2,0) in the flavor and color 
spaces. 

Introducing the auxiliary fields 


with i(Dp“'(x) = iy^b^ — m + fisyo + Bsyoys and the second 
Gell-mann matrix A 2 in color space. After integrating out the 
bispinors and eliminating the double counting, the thermody¬ 
namic potential Q = -T/V\nZ can be expressed as 


cr - -2Gsipp, 
n - -AGsipiysTip, 

Ac = - 2G Blip cfseoysip , 

A* - -AGoiipsesysipc, (23) 


the Lagrangian density becomes 


£ - ip[irb. 


mo- cr - iysT ■ n + jigyo + Bsyoys \ P (24) 


-- [Acipsesyspc + A^ipcECoysp] - 


4Gd 


On — mo)^ Ac i 


4 G 5 4G/) 2 V 4 


For q - 0,1, the trace can be worked out easily 


'.;=± 


+ ^ +T\n[l ( 31 ) 


n,P3 


In mean held approximation, supposing constant condensates where 'Zjn,p, - ff 2n=o f ^ with a,, = (2 - b„o)/2 is the 
(cr) - m - mo, (tt) - 0 and (Ac) = (A^ = Ac, we have the summation over the Landau level and the integration over the 
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longitudinal momentum in the external magnetic field, and the 
quasiparticle energies are defined as 

E^p) = + ipsf + m?- , 

Elipi) = + ipsj +iv?-. (32) 

For q - ±112, the trace can be evaluated exactly, due to the 
same sign of the psyajs terms and the same electric charge in 
the diagonal terms of 5'^. Following a similar procedure as in 
Ref. IS, we get 

0 +1/2 + 0_i/2 = - Z Z K' + 2F In (l + (33) 

/,y=± «,P3 


with the quasiparticle energies 


EUp^) = 



(34) 


For simplicity, we consider in the following only the case 
at zero temperature T - Q and in chiral limit with mo - 0. In 
this case, the thermodynamic potential can be simplified as 


Q. 


E? 


f 


d^p 


Eb, 


■Z 


E" 


^-e{E^;-PB) + Y^PB-E'f) 


-^e{E'i-pB) + ^e{pB-E”)±Ei 



with the unit step function 6{x). To remove the divergence 
in the integrations, we use here the Pauli-Villars regulariza¬ 
tion scheme fiH, which is much softer than the hard three- 
momentum cutoff and works better in the case with a nonzero 
chiral imbalance, in comparison with the vacuum regulariza¬ 
tion used in Section|III] From the numerical results shown be¬ 
low, the Pauli-Villars regularization will greatly suppress the 
artificial oscillation of the quark mass in external magnetic 
field d. In this way, the regularized version of the thermo¬ 
dynamic potential is given as 


Q = 


AGs AGd 


i.i=± 1=0 


f 


d^p 

(2;r)3 


%(8ypB) + Y(^(PB-&;) 
Pb . 


•z 




(36) 


with fi - 31^-61-1-1 and the regularized quasiparticle energies 

£°(l,p) = ^j{E°{p)f + 1 A 2 , 

SliUps) = ^j{E'/(P3)f + 

&ljil,P3) = ^{&f{l,P3) + jpBf + ^l (37) 


The chemical potential and magnetic field dependence of 
the order parameters m{pB,ps,B) and Ac{pB,P 5 ,B) is con¬ 
trolled by the minimum of the thermodynamic potential 
dill dm - 0 and diljdAc - 0 which lead to the two gap equa¬ 
tions 


m 

2 ^ 


A, 

2Gd 


•ZZf'lf 

ij=± 1=0 ^ 


d^p 1 
(2;r)3 2£| 

S n/A., 

i +JBb 


oOiSi-EB) 

■njl 


sys'ij 


^ZZZ-ir: 

ij=± 1=0 n,p2 iJ 


(38) 


Note that with the Pauli-Villars regularization we used, the 
thermodynamic potential (l3^ is still divergent, but the gap 
equations ( l38l l we are interested in are convergent. In fact, 
the quantity Q. corresponds to the pressure except for a sign, 
and only the pressure relative to the vacuum can be mea¬ 
sured. While Q.{pB,P 5 ,B) is divergent, the physical potential 
Q.(Pb,P 5 , B) — Q(0,0,0) is convergent. 

It is easy to see that a trivial solution of the gap equations is 
m - Ac - 0, which corresponds to the phase with both chiral 
and color symmetries. The phase with both chiral and color 
symmetry breaking described by m,Ac A 0 does not exist, 
due to the conflict between the two gap equations in chiral 
limit. The possible nontrivial solutions include the normal 
phase with m 0, A^ = 0 and the color superconductor with 
m = 0, Ac 0. 

We choose the Pauli-Villars parameters as A = 0.859 GeV, 
GsA^ - 2.84 and Go - ^Gs IA idUli^l . Let us first look at 
the chiral imbalance effect on color superconductivity, shown 
in Fig|7]at pB — 0.4 GeV and B - 0. Since the diquark con¬ 
densate can be expressed in terms of chiral quarks as 


Ac = -2GDi{<Pc>se3j$>l/) 

= 2Gd [<iAlr2roee3r5iAL> + <iARr2roef3r5<Afi>] > (39) 

the two condensed quarks are with the same chirality, and 
therefore the chiral imbalance should enhance the color con¬ 
densate. 

The baryon chemical potential dependence of the two order 
parameters m and Ac of chiral restoration and color supercon¬ 
ductivity at fixed chiral imbalance and magnetic field is shown 
in Fig0 The system is in normal phase with chiral symme¬ 
try breaking {m + 0) at low pB and color superconductivity 
phase (Ac + 0) at high pB- The chiral condensate keeps as a 
constant in the normal phase, suddenly jumps down to zero at 
the critical baryon chemical potential and remains zero in 
the color superconductivity phase. In contrast, the color con¬ 
densate remains zero in the normal phase, suddenly jumps up 
to a nonzero value at and then increases with //b in the 
color superconductivity phase. The two phase transitions at 
the critical point p'^ are both of first order. With increasing 
chiral imbalance, both the chiral and color condensates are 
enhanced, but the critical point shifts towards right. This in¬ 
dicates a stronger catalysis effect of chiral imbalance on the 
chiral condensate than on the color condensate. 
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FIG. 7: The color condensate A^. as a function of chiral imbalance 
/J 5 at fixed baryon chemical potential fig = 0.4 GeV and vanishing 
magnetic field B = 0. 



FIG. 8 : (color online) The chiral and color condensates m and A^. 
as functions of baryon chemical potential yUg at fixed magnetic field 
eB = 0.3 GeV and chiral imbalance ps = 0 (solid lines), 0.2 GeV 
(dashed lines) and 0.4 GeV (dot-dashed lines). 

The magnetic field effect is shown in Fig|9] While the 
de Haas-van Alphan oscillation for quark mass m is not ob¬ 
served, which is consistent with the results in Refs, llz^ 1^ 
I 43 I] . the oscillation shows up for color condensate at any 
combination of chiral and baryon chemical potentials, similar 
to the results in Refs. iflOl - U^ . We would like to point out 
the similarity between the two chemical potentials in the gap 
equation for color superconductivity. In this phase with m - 0 
and Af 0, the excitation energy becomes 



and we can introduce two combined chemical potentials = 
\fi 5 + fisl instead of fiB and In this way, there should be 
no difference between the case with (jj.b + 0 ,fi 5 = 0) and the 
case with {fiB - 0 ,/j .5 + 0). The discrepancy between the 
line with {^b,Bs) - (0.4,0) GeV and the line with {fiB,B 5 ) - 
(0,0.4) GeV shown in the lower panel of Fig|9]is due to the 
Pauli-Villars regularization which makes jj .5 different from 
by introducing two large masses A and V2A. However, the 
regularization mainly leads to a global shift between the two 
lines, it does not change the structure of the magnetic field 



FIG. 9: (color online) The chiral and color condensates m (upper 
panel) and A^. (lower panel) as functions of magnetic field B at fixed 
baryon and chiral chemical potentials pg and ^ 5 . 



(GeV) 


FIG. 10: (color online) The phase diagram of chiral restoration and 
color superconductivity in pg - S plane at fixed chiral chemical po¬ 
tential /J 5 . 


dependence. 

Finally, we show in Fig[T0] the phase diagram of chiral 
restoration and color superconductivity in the plane of baryon 
chemical potential and magnetic field at fixed chiral imbal¬ 
ance. As can be seen, the de Haas-van Alphan oscillation for 
the critical baryon chemical potential shows up clearly at 
vanishing chiral imbalance. Such oscillation for chiral phase 
transition has been observed and well explained in Ref. i^. 
At large chiral chemical potential - 0.4 GeV, the critical 
baryon chemical potential is almost a constant = 0.35 GeV, 
the oscillation is washed away by the strong chiral imbalance. 
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V. CONCLUSIONS 

The phase transitions with chiral imbalance and magnetic 
field in strongly interacting quark matter are investigated at 
finite isospin and baryon densities. In the frame of extended 
NIL models, we focused on the magnetic effects on the pion 
superfluidity and color superconductivity and the related de 
Haas-van Alphan oscillations. For pion superfluidity, due to 
the problem of mixed quark propagators at different Landau 
levels, we treated the chiral imbalance and external magnetic 
field separately and take into account the interaction between 
the charged pion condensate and the magnetic field in the GL 
approach. For color superconductivity, we self-consistently 
studied the chiral imbalance and magnetic field effects on the 
order parameters of chiral restoration and color superconduc¬ 


tivity. 

For the pion superfluidity, the critical value of isospin 
chemical potential is suppressed by the chiral imbalance but 
enhanced by the external magnetic field, indicating respec¬ 
tively a catalysis and an inverse catalysis effect. The latter 
is consistent with the lattice simulated magnetic effect. For 
the color superconductivity, the chiral imbalance leads to a 
weaker catalysis effect compared to chiral symmetry break¬ 
ing, and the magnetic field results in a de Flaas-van Alphan 
oscillation on the phase transition line. However, the oscil¬ 
lation is washed away when the chiral imbalance is strong 
enough. 
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